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Abstract 



We pointed out in this work that in dealing with the BEC-BCS crossover 
problem, the dynamic effect of <b(t)>is not negligible. Accordingly, an equation 
of motion approach was devised to calculate the Green's functions. Based on 
our result, we concluded that instead of crossover, BCS states and Bose- Einstein 
condensation always coexist. 
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The physics of BCS states in a dilute system was first studied by Eagles [1] 
who proposed the possibility of forming Cooper pairs before the onset of super- 
conductivity. The BEC-BCS crossover problem was discussed by Leggett[2] in 
1980. He pointed out that there is no distinct boundary between Cooper pairs 
and Bose-Einstein condensation. As elaborated later by Nozieres and Schmitt- 
Rink[3], in the strong coupling or low density limit, fermions will combine to 
form bosonic molecules and condense at sufficiently low temperature. As the 
coupling is weakened or density increases, the molecular wave functions over- 
lap. With sufficient overlapping, the molecular BEC will be transformed into 
BCS states. The authors [2,3] saw no phase transition. In their opinion, the 
amplitudes of the states of phase coherence (BEC and BCS) undergo smooth 
variation. 

On the other hand, there has been marvellous experimental progress[4-8] 
recently. Fermionic atoms in a trap can form bosonic molecules via Feshbach 
resonance [9]. The coupling strength can be adjusted by changing the applied 
magnetic field. The energy gap was found on both sides of resonance and its 
spectrum is compatible with theoretical calculation[10]. Hence, the scenario en- 
visaged by Leggett[2] and Nozieres and Schmitt-Rink[3] seems to be realized, 
with one complication: the fermions and bosons have internal degrees of free- 
dom. They are the spins of electrons and nuclei. The above mentioned Feshbach 
resonance is fulfilled with the aid of hyperfine interaction. As a result, the state 
of electron spins of the atoms in the open channel are different, in general, from 
those in the closed channel. 

Much theoretical work has been devoted to this subject. Kokkelmans et. 
al.[ll] made a thorough analysis of Feshbach resonance of the cold atom system. 
Their deduction greatly simplified further computation. Ohashi and Griffin[12] 
adopted the orthodox Feynman diagram approach. Stajic et' al.[13] applied the 
concept of pseudogap originated from high Tc superconductivity to explain the 
experimental results. In this paper we would like to raise a neglected point 
which has profound effect on one's understanding of the BEC-BCS system. 

We started with the of fermion-boson Hamiltonian[14]: 



where c Pi<T and 6 q are the fermion and boson operators, s p and w q are their 
energies, v is the detuned energy of bosons, and a is the spin index. The terms 
in the second line denote the combination of atoms and dissociation of molecules. 
They are necessary because, as we stated before, the states of spins had been 
changed. Cooper pairs have different spin states from those in BEC. So they 
should be treated as different species. 

Now we calculate the Green's functions G a (T,T l ;k) = — (Tck, CT (r)cj c a (r')) 

and F*^(r,T';k) =(Tcj c Q (r)cj c ^(r')) where T is the time-order operator. The 
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approach of equations of motion was applied. 

J-G T (T,T';k) = - S (t-t') -s p G^r,r';k) - [C/F iT (r, r) - g(b )]F^(T, r'; k) 

+ .g^(T[6 q (r) - (6 >]c t _ k+q ,i W4, T (^)> (2) 
q 

and 

(t, r'; k) = £ P -F|*j (r, r'; k) + (r, r) - 3 (fe )]G i (r, r'; k) 

+ g]T<T[6t( T ) _ (6t)] Ck+qii ( r ) c t )i(T ')) ( 3 ) 

q,er 

where F a p(T,r) = X^ P ( c -p,a( T ) c p,/3( T ))- Here we assumed the Hartree-Fock 
terms can be absorbed into the single-particle energy. A term (& q =o) had been 
intentionally added to and substracted from the right hand side to facilitate 
later decoupling. 

While developing superconductivity theory, F a p{T,r) is usually set to be a 
real constant independent of r. The prescription is to subtract N/j, from [15]. 
Similarly, in calculation involving BEC, (6o(t)) is also treated as a constant. In 
the case when both F a p(r, r) and (&o(t)} are present, this method in general, 
can only make either one of them a constant but not both. This point was 
ignored by the preceding works on the BEC-BCS crossover problem. In this 
paper, we chose to let F a p(r, r) = F a p be a constant and (bo(r)) a function of 
r. This is the gist of this work. 

With H replaced by H — N^i 7 we have to substitute £ k = £ k — ^ for e k and 
v' = v — 2/i for v in cqs. (1-3). If (&o(t)) was treated dynamically, we arrived 
at the next set of equations of motion: 

^<T[^(r) - <6o>]cL k+q>i (r)4 >T (r / )> 
= [-K + ,')H N ](T[Mr) - (&o)]c t _ k+q!i (r) C t T (r')) - K + ^)(&o)(TcL k+q;i (r)c k;T (r')) 

+ £<T[6 q (r) - (foo)]^ +q ^ T (r) C t _ k+q _ q , a (r) Cp , T (r) Ck!T (r')) 
p,q' 
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and 

|:(T[6t( r )-(6t )]ck+q;i(r)c t a(r0) 
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+ f/^(f [6 q (r) - (4)]c P+q ,, T (T)c k+q , +q;i (T)c p , T (T)4 ii (r')) 
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At this stage we applied the decoupling between the fermion operators and 
boson operators. As before, we assumed the contribution of the interaction 
terms of eqs. (4) and (5) (the third and fourth lines) can be incorporated into 
the one-particle energies £ k and v. The resulting forms will be comprehensible 
and easy to analyze. Thus, eqs. (4) and (5) became 

[^+(w q +i/)-&- q ]<T[6 q (T)^ 

(6) 

and 

[|--K+^+& +q ]<T[& q (T)-<&t)^ 

(7) 

Taking Fourier transform of eqs. (2), (3), (6) and (7) and making substitutions, 
we obtained 

G k , T «w„) = — — , (8 

and 
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(11) 

where B = + (6c - ^') 2 + (^fl ~ S^) 2 ]/ 2 - Clearly, the Green's functions still 
possess the conventional form of BCS theory. However, we accounted for the 
most important ingredient: the coupling between F*^ and <f), to insure that we 
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got a valid physical picture. In fact, as shown below, our results are qualitatively 
similar to those of ref. 12 who had considered many Feynman diagrams. 

In order to study the resulting physical properties, we followed Eliashbcrg[16] 
to write down the form of the thermodynamic potential 

« = ]T£k+(L/F iT -. 9 ^ 

k n.k q 

(12) 

Here, we have kept only the zeroth and first order terms. The relation F*^ — 
is also helpful. Now we take the variational approach. 

N = -™ = V[l + ^ tanh^ + ^tanh^i - tanh^^] 

+ 20 2 + 2^Tn B K), (13) 
q 

and 

— = 24>v - gF u - iJ-Qj- tanh — + — tanh —] = (14) 
where ngK) = l/[e^ Wq+l/ ) — 1]. Additionally, we have the gap equation 

Ftl = h U = \ 2E + ^-^) 

UF^[E*_ M - (& - v'f] - g^[El M - & (& - v')] 
V 2f?_, k (^ ik -^ ik ) 

(15) 

With eqs. (13-15) we can solve for three unknowns: fi, <j) and F*^. 

We showed the results in Figs. 1 by plotting <f> and FZ, versus temperature 
at detuned energies v = ±0.5Ep in the solid and dashed lines respectively. 
The parameters are the same as those in ref. 12, i. e., g = —0.6Ep, U = 
0.3Ep and a cutoff factor exp[— (e^/2Ep) 2 ]. In Fig. 2 the amplitudes of the 
states of phase coherence F*^ and 4> at ksT — Ep/3 (in the dashed lines) and 
ksT = Ep/6 (in the solid lines) are plotted against the detuned energy. It is 
remarkable that in Figs. 1 and 2 FZ, and <j> coexist in the whole range and vanish 
together. There is no crossover. In fact, similar results had been obtained by 
previous calculations[ll,12] with the static approximation (b(r)) ~ <p and the 
resultant relation <j> = —gFZJ{v — 2p) and energy gap A = [U +g 2 /{v — 2fi)]FZ^. 
Note that one of the consequences of the approximation is that <j> and F*^ are 
proportional. However, our approach did not have the proportional relation a 
priori. Therefore, we are able to conclude that the original concept of BEC-BCS 
crossover cannot be applied here. 

In Fig. 3 we presented - ImF^(B) versus EjEp. On the BCS side [y > 0), 
the feature of Cooper pair density of states is clearly shown. It gradually became 
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molecule- like when v became negative. The reason, as pointed out by Leggett[2], 
is that when the chemical potential becomes negative, the divergence in the 
spectrum disappear. Note that in the whole range, the amplitude of BEC is 
greater than that of BCS due to the parameters we chose. Yet the spectrum 
underwent drastic change. Clearly, the shape of spectrum is not an indication 
of whether it is BEC or BCS state. 

It should not be surprising that BEC persists deeply into the BCS side of 
Feshbach resonance. Although the detuned energy is positive, forming molecules 
can still be energetically favorable. It is because most atoms (fermions) have fi- 
nite amount of kinetic energy. On the other hand, on the BEC side, it seems that 
molecule formation is advantageous. However, BCS state is helped by having 
the energy gap. The larger the gap, the lower the energy. For negative detuned 
energy, the chemical energy is also negative. As pointed out by Leggett[2], the 
gap under this situation is at least y/ fi 2 + A 2 . The energy gain by forming 
Cooper pairs grows with the magnitude of the chemical potential and the de- 
tuned energy and thus, the coexistence. Furthermore, if we simplified our model 
by considering only the condensed bosons (6 q= o), then the Hamiltonian in cq. 
(1) would remind one of the Anderson model[17] or Fano resonance [18]. The 
eigenstate should be a mixture of localized state (condensed bosons) and con- 
tinuum (fermions). This analogy, though not completely compatible, indicates 
that the states of phase coherence should have two components. In conclusion, 
we pointed out that cither (6(r)) or F*^ (t, t) should be time-dependent and 
our calculation show that BCS state and BEC always coexist. 
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Figure captions 

Fig. 1 F*^ and <j> versus temperature at v = 0.5Ep (solid lines) and v = 
—0.5Ep (dashed lines). 

Fig. 2 and cf) versus the detuned energy v at temperature at Ep/6 (solid 
lines) and Ep/3 (dashed lines). 

Fig. 3 — ImF^(E), the spectrum funtion of Cooper pairs versus E/Ep at 
different detuned energies. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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